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Abstract 

A structural expansion for the static ground state energy 
of a simple metal is derived. Two methods are presented, one 
an approach based on single particle band structure which 
treats the electron gas as a non-linear dielectric, the other 
a more general many particle analysis using finite temperature 
perturbation theory. The two methods are compared, and it is 
shown in detail how band- s true ture effects, Fermi surface 
distortions, and chemical potential shifts affect the total 
energy. These are of especial interest in corrections to the 
total energy beyond third order in the electron ion inter- 
action and hence to systems where differences in energies for 
various crystal structures are exceptionally small. Preliminary 
calculations using these methods for the zero temperature 
thermodynamic functions of atomic hydrogen are reported. 


I. INTRODUCTION 


Recent work in the theory of metallic phase stability 

has met with moderate success in accounting for the most stable 

crystalline structure, binding energy, and compressibility 
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of a simple metal. * The theory depends upon a perturbation 

expansion of the ground state energy (T = 0°K), usually to 

second order in the Fourier components of the pseudopotential 

evaluated at reciprocal lattice vectors. In certain cases, 

however, the energy difference between structures is so small 

that it is essential to consider higher order terms in a 

structural expansion for the energy. A case in point is 

atomic metallic hydrogen for which a second order calculation 

of the ground state energy per proton using a RPA dielectric 

function gives (static lattice) energies of -1.01532, -1.01597, and 

-1.01537 Ry respectively for the SC, FCC and BCC structures at a 

density (r =1.6) near the zero pressure metastable equilibrium, 
s 

The procedures for constructing the perturbation expansion 

3 

have been known since 1958 when Hubbard developed a diagrammatic 
technique based upon solutions of a one-electron Hartree-like 
equation, a method which ultimately enabled him to express 
the energy in terms of the solutions to an integral equation. 

Later, self-consistent methods were proposed by Cohen^ who 
treated the ground state properties of a solid along the lines 
of the dielectric formulation of Nozieres and Pines'* for the 
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electron gas. More recently Brovman et al.^ have used a 

modification of Hubbard's technique to calculate both binding 

energies and phonon spectra for simple metals. Lloyd and Sholl^ 

have also presented explicit expressions for third order 

corrections to the total energy using an analysis similar to 
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that of Hohenberg and Kohn , and Harrison has discussed the 
interpretation of these contributions in terms of three body 
interactions. What we present here is an explicit structural 
expansion which is convenient for calculation of ground state 
energy as a function of density and which is simply related 
to the eigenvalues of the one-electron band Hamiltonian. We 
shall discuss its relation to a more complete solution given 
in terms of the T = 0°K limit of finite temperature perturbation 
theory. Finally we shall discuss certain differences between 
the present work and the previous theories mentioned above 
and apply these techniques to a calculation of the ground 
state properties of atomic hydrogen. A comprehensive Bravais 
lattice survey of the binding energy to third order in electron 
ion interaction for this solid has been carried out by Brovman 
et al.^ . The purpose of our calculations is rather to study 
the magnitudes of higher order corrections, in support of which 
we shall present numerical values for SC, FCC, and BCC lattices. 


II. FORMULATION OF THE PROBLEM 


We consider in this section the problem of computing the 
total energy of a system of N interacting electrons in a 
static periodic one body potential. Later^ we shall relax 
this restriction and consider the modifications arising from 
dynamic effects. To begin with we shall restrict our con- 
siderations to T = 0°K and subsequently extend the analysis 
to non-zero temperatures. 

The Hamiltonian for our system thus restricted may be 
written as 


H = H + H . + H. . 
ee ex 11 


( 1 ) 


where H describes the kinetic and interaction energy of a 
66 


system of coupled electrons, i.e 
H ee = S l P i + 4 l 




* 2 

|r.-r . \ * 

'~i ~J 1 


( 2 ) 


H^ describes the interaction energy of the rigid lattice of 


ions of valence Z, i.e. 


H. . 

li 



(3) 


and H . describes the interaction of electrons with the lattice, 
ei 

i.e. 


H . 


ei 



( 4 ) 


In (3) W(R,,R 0 ) is the bare ion-ion interaction, and in (4) 

^4 
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V(r) is the periodic one body potential. We may express H 

rw« 

in terms of second quantized operators, i.e. 

1 V 2 h 2 c.2 + 

2m l P i 2m 1 C k C k 


x 

i 


*1 |r.-r.l -*1 ’ 

Ij k.k’. S 


c k+q C k-q c k ,c k 


ee 


(5) 


( 6 ) 


where 


ftw(q) = J e 
ft 


2 2 
-iq*r e ,3 4ne 
n — d r = — 7T- t 
r 2 

q 


(7) 


is the Fourier transform of the bare Coulomb interaction (ft 

12 

being the volume of the system ). There is the usual problem 
of handling the q = 0 term. To resolve it we carry out the 
following series of manipulations (the thermodynamic limit 
being taken as the ultimate step). First, we subtract from H 
the term 

^e 2 (f) 2 JIA , (8) 

ft 

that is, the interaction term in H is replaced by 

t 


ee 




*1 

1 „1 ' 


1 .3. 


(9) 


ft 


and H accordingly becomes 
66 


H e e “ 2S 1 h2k2 Vk + * I w ( S>4+/k' -q C k' c k 

A/ ^ rv/ r3 f*j 


( 10 ) 


k,k f 


which is the familiar electron gas Hamiltonian, denoted in what 
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follows by H . We now add (8) to H. . . Thus 
-'eg 11 

' P I 2 . o 

H ii-I “<W + 4 J d r l d r 2 ’ 

ap Cl 


(ID 


where P = — . The term which has been added is the self- 
o i 2 

energy of a uniform background of negative (or positive) 

i 

charge. To (11) we add the interaction energy of the ions 
with this negative background so that H. . becomes 

■ p 2 

H ii - * l w <W + i It: d3r i d3r 2 - 1 1 p o w( £>V d3r > (12) 


which, if we neglect Born-Mayer terms, is just the Madelung 
energy for the assembly of ions. Finally, we subtract the 
same interaction energy from the last term in H, H^, obtaining 


:i -I v <Ei> + IJ p o w( E>5d )d3r • 


(13) 


The original Hamiltonian has now been separated into three 
well defined parts. Taking its average over the ground state 
we have as the expression for the total energy per electron. 


I = H { < H eg > + < H ei >} + =M (14) 

where is the Madelung energy, i.e. the energy per electron of 
a lattice of positive ions in a uniform background of negative 
charge. Notice that the first term in (14) is not the energy 
per electron of the interacting electron gas since the ground 
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state wavefunction is that appropriate to an electron gas in 
which a periodic array of ions is immersed. 

Let us consider the second term in more detail. For 
even a simple metal the interaction potential V(r) is not 
known in general from first principles. From the point of 
view of band theory, however, it may be well represented by 
a weak pseudopotential, at least for the valence states. 

(We set aside in this discussion questions of core level shifts 
and their effect on the total energy.) If we make this pseudo- 
potential approximation and furthermore consider a local 
approximation in which the periodic potential is a simple 
superposition of bare pseudopotentials at each lattice site, 
then (13) becomes 


ZeP 


H . = Y v(r.-R ) + Y [ -j — r d r 

ex L ~x ~a L J r-R 

. „ ~a 

ia a ft 

or in terms of the Fourier transform of v, 


H . 
ei 




ik. r-? 
e ~ 


a ft 


ZeP 
o 

r 



(15) 


(16) 


where 


ft 


v(k) = f v(r)e d^r , 

r+s J r«-» 


ft 


and 


P< i} - X . 


(17) 


a 


(18) 
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In particular, the k = 0 term is given by 


Lim 
k 


0 [NiNv(k) + N.ZeP o J ^ ] * 


( 19 ) 


Q 


where fL is the number of ions, N = ZN^. As an example, for 

a potential which is Coulombic beyond a certain "core*' radius 
13 

r and vanishes otherwise we have 
c 

o r 

60 .3 _ c 


„ Lim .. N 2 I* d r p f * 

a k-o v( ~ } “ - Ze I — + J v( £ ) 

~ r o 


d 3 r 


( 20 ) 


Hence the long range parts in (19) cancel and we are left with 


N, V c 


H ei = l p k v(k)e 1 ~'£i + N ■£ J v(r)d 3 r + N2eP o J 

MO ~ 0 c 

1 

which we rewrite as 

H . = V P^ v(k)e 1 ~*~i + NE , 
ei L k ~ c 

• i ~ 

x,k 


c ,3 


dr 


( 21 ) 


( 22 ) 


14 


where the "core" contribution E c is independent of structure, 
and the prime means that k = 0 is excluded from the summation. 
Thus the ground state energy (T = 0°K) can be written in 
the form 


I 

I " S < H eg » + I < l P k 1)v( !S )P -k > + E c + > 


(23) 


o xk • ^ « 

where P_^ = ^ e 1 • For a lattice of bare protons, we note 

^ • 

that (23) is 1 exact with E = 0 and v(k) = w(k). However, for 

C 0SJ 


the general case it is approximate since it is not clear that 
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a single-particle equation describing the band structure with 
a local v(r) can be derived from H as given in (23) with the 


same v(r). Moreover, it is not strictly correct to write the 
lattice potential as a simple superposition. With these 
reservations, we may address ourselves to the task of computing 
the average 


- < h + y py^voop , > 

N eg L> k ~ -k 


(24) 


n f X ) 

If we treat = ) P^ v(k)P ^ as a perturbation, then 

. r*J 

k 


the unperturbed problem is. the interacting electron gas. 
Indeed, the problem is that of a dense distribution of identical 
impurities in the electron gas except that for a crystal, the 
impurities are arrayed in a definite order. Alternatively, 
one may simultaneously treat both electron-electron and electron- 
ion interactions as perturbations and carry out the usual 
double perturbation expansion. In the following sections we 
present two methods for computing the energy shift due to H^, 
one closely related to a single particle picture, the other a 
more general many particle method. 



III. BAND APPROACH 


In this section we consider the calculation of the ground 
state energy from a single particle point of view. The physical 
picture is the following. We have a system of electrons whose 
interactions with the static lattice are described by a pseudo- 
potential. The electron gas may be viewed as a non-linear 
dielectric and the pseudo- ions as the source of external 
potentials which induce charge density responses in it. The 
energy associated with this induction process is given by the 
well known expression"^ 

6W- j6V(r)P(r)d 3 r (25) 

This is the work which an external contrivance must do in 
changing the potentials from some value V to V + 6V. In 
terms of Fourier transformed quantities this becomes 

6W = SI Y 6V(-k)P(k) . (26) 

/ . o/ 

k 

The contribution of the electron-ion interaction to the total 

energy is then given by 
V 

W = J 6W (27) . 

o 


In general, we may write the charge density P(k) as 
P(k) = X 1 (k)V(k) +X X 2<B>3) V <S+3>VC- 3 ) 

3 < 28 > 

+ I x 3 ( b3i’32 )v( h + 3i + E2 )v( "Si )v( '32 ) + •••> 

Si >32 
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the first term of which is the usual linear response expression. 
In appendix A we show that this leads to an expression for 
the change in energy given by 


n ^ ■ w " *1 Xi(k) v (-k> v (k) + x 2 ( ~’3 )v( !3 + 3 )v( "3 )v( ‘y 

k k,q 

rv rO 

+ * l x 3 ( S>3i-32) v ft + 3i + 32 )v( -3i) v <-S2> v (-y + •• 
b’3i>3 2 


(29) 


which we shall refer to as the band structure energy and 
which is determined from the induced charge density through 
(28). Note that V(0) is to be excluded from the summation 

rv< 

(a requirement of charge neutrality as discussed in section II). 
Equation (29) thus presents us with a well defined method for 
calculating in terms of the charge density. 

From the point of view of single particle band theory we 
calculate the charge density from the Bloch wave function 
of an electron in a periodic potential V. In terms of plane 
waves 

< r|k-K. > = ft ^ (30) 

»s/ 

the wave function is written as (we assume a Bravais lattice) 

“I c k-K s I c il i:> > 

~ 17 ~ ~ 


where the coefficients c^ satisfy the equations 


(31) 
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(£ - E + V ) c + V m c- + y .c. = 0 

v o oo o 01 1 L oi l 

i^O, 1 

^10 c o + «? 1 - E + ^l )c l + l V li C i = 0 

i^O, 1 


( 32 ) 


c + + (cf . 

xo o xl 1 X 


E + V. .)c. + Y V. .c. 
xx' x L xj j 

j^0,l,i 


= 0 


wx 


- Z. r) 

th £ . = — (k-K.) and V. . = < k-K. I'vlk-K. > . An iterative 
x 2m xj x 1 


' 1 ' xi ~ ~x 

solution of equations (32) yields a Brillouin-Wigner expansion 


for the namely: 

C i = c o{e^T7 + I 




V..V. 

1-1- -1.0- 


(E- <f . ) (E- £. ) 


>1 

jk 


-ij . 1 k . jSP., 


(E- <5 i )(E- cf. ) (E- <f k ) 


+ . . . 


’ V. 

i-1 il + V 11 ik kx 
\) 2, (E- C 


(33) 


, (E- cf . ) (E- cf 

1 j 1 3 jk 


.)(E- (f j ) (E- cf k 


) + ...} 


where the prime excludes 0,1. Equations (33) lead to a 
folded secular equation 

(cf - E + U qo )c o + U Q1 c 1 = 0 

U 10 C o + ^ 1 " E + U 11^ C 1 = ° 
with the U's defined by 

/■O * f>- 


(34) 


V..V. 
Ax xm 


' V..V..V. 

Ax ij im 


(35) 


it = t; i V Ai im , V 1 i i ... | 

Am V Am L E - <f . L (E- . ) (E- £ . ) 

i 1 ij 1 3 

In (35) the prime excludes A, m from the summations. Note 

that although U„ = U U„ n ^ U„ . The folding transformation 
& Am mA A-m, 0 Am b 


-13- 


is valid for any l, m and accordingly, 

(<£ 4 + U.„ - E)c . + U„ c =0 
i Hi H im m 

U c + (cf + U - E)c =0 
mj l & m mm m 

These equations define a two band (upper and lower ^ ^ ) 
situation for which the solution for K, ^ 0 is: 

E, x = <5 + u ( - } + |U ( “ ) |{y ( “ ) - (1 + (v ( - ) ) 2 )H 

(-) o oo 1 om 1 l m v v 'm J 


(36) 


(-) „ 


m 


m 


(-) _ 


7-T- { 6 - £ + - u (_) } 

(-) , l m o mm oo J 

Of + u (_) 


aju''"' 1 

1 om 1 


u<-V-> 

mo o 


(<f + U (_) - E, O 

m mm (-) 


U 


(-) 

om 


(37) 


hi)l c (-> 

o 


or 


(“) 


m 


+ (v<->> 2 ) 4 } c<-> 

om 


(38) 


A similar expression holds for the upper band with (-) - (+) 

atld { Y m _) " (1 + { Y m +) + (1 + (Y m +))2) ^} * We ma V 

use these results to calculate a charge density, i.e. 


*£> = 2 1 I c i c i < ilj >< £lj > > ' 

[k] ij 

where ^ denotes a summation restricted to occupied levels, 

[k] 

The Fourier transform of (39) gives 

p *=n ivw 

fk] i 


(39) 


(40) 
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which for a single occupied band reads 

p 2 v v 

^ i, t I {v<^ 2 > 4 } 2 } ' 

i¥o 

Alternatively, this may be rewritten using (37) as 


( 41 ) 



r 2U QJ& , y U i-l,o U i,o 1 

*<V <-W ±%.i (E r £ i-r u i-i,i-*><V 


{'+1 


|U, 


12 - 


i*o <Y W 


} 


(42) 


These last two expressions are easily generalized if two bands 
are occupied. If more than two bands are occupied it is 
necessary to begin with the folded secular equation appropriate 
to that number. We note again that in (41) the k summation is 
only over occupied levels. Thus we are summing up to the 
true Fermi surface rather than within a Fermi sphere (the 
more common situation in perturbation theory). 

The above expressions although formally exact within the 
one electron approximation are difficult to use in practice. 

If we knew the analytic dependence of the U s on V we could 
perform the integration in (27) (for example, by associating 
with V a coupling constant over which we ultimately integrate) 
and finally carry out the sum on k. However, only in the 
extreme approximation of retaining a single y is this analytically 
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tractable. We can, on the other hand, expand the expression 
for P in powers of V. If we then assume that V(k) = V(k)/e(k) 

Xj ~ ~ ~ 

where e (k) is the static limit of the electron gas dielectric 
function, it is possible to calculate the energy shift from 
(29). The results of such an expansion are given in appendix B. 
In the following section we shall derive an expression for 
the energy shift from a more complete theory and see that the 
simple theory above must be only slightly modified. 


IV. FINITE TEMPERATURE RESULTS 


In this section we calculate the total energy of the system 

of electrons and ions using the techniques of finite temperature 

perturbation theory. If we choose as the unperturbed system 

one having a spherical Fermi surface (e.g. non-interacting 

or interacting electron gas), it is, in fact, necessary to 

use this method, a consequence of the fact that for interacting 

electrons in a periodic potential, the adiabatically generated 

state of the zero temperature method is not the true ground 

state, no matter how weak the lattice potential. The 

state generated adiabatically from a spherical ground state 

can- -never depart from a state with a spherical Fermi surface 

17 18 

and cannot produce the crossing of levels ’ resulting from 
the imposition of a periodic potential. In the finite tempera- 
ture theory, however, the mean occupation number of a given 
quantum state is no longer restricted to be either 0 or 1. 

Thus the Fermi surface of the unperturbed system is permitted 
to distort in such a way that the thermodynamic potential is 
minimized subject to the constraint of fixed overall density 
and in consequence the true Fermi surface is attained at 
each stage of the calculation. 

The temperature formalism is most simply stated in terms 

of Green functions. We shall follow the exposition of Martin 
19 

and Schwinger and define the single particle Green function as 
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where the brackets denote the grand canonical ensemble average 


< <3 > = 


Tr 


Tr e 


-p (H-uN) 


(44) 


are Heisenberg field operators, a, p are spin indices and 
T is the time ordering operator for real t (and the 'it' ordering 
operator for imaginary times.) We Fourier transform G and 
write the result itself as a Fourier series 

G ae ( £l> E 2 ;t) S i Jd 3r L d 3r 2 e'^rEl e ^rZ2 G (r^r^t) 


a 


n 


G a S ( El>E2 ;t) 


(45) 

(46) 


V 


where a> - 


TT 


(2v + 1) + p , v = 0, +1, . . . , so that 


v (-ip) 

(-ip) 

G aB ( El’E 2 ;“v ) = J G aB ( El’E 2 5t)d,: • 


(47) 


These results are consequences of the boundary condition 
satisfied by G for imaginary times. The average value of a 
one-body operator is given in terms of the Green function by 

< v > = “ Y V(-p)ftG ; (k,k-p;cp )e^° , (48) 

P Z_J OL (X r*> r*J \J 

jg»k> v 


(where for the sake of brevity G = ) G ) . 

J aa L aa 
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In order to compute the ground state energy we use the 

20 

statistical mechanical theorem which states that for any 
parameter X in the Hamiltonian, 


M = < > 

ax ax 


( 49 ) 


where S is the thermodynamic potential, the differentiation is 
at fixed T, V, p, and the average is that defined in (44). 

For the Hamiltonian we take that given in (24). If we associate 
a coupling constant X with the bare interaction V, we then have 
upon integration. 


SGa) = 8 0.) + f f < XV • 

o 

To calculate the ground state energy we take the T = 0 
limit of E(p) + pN, i.e. 


(50) 




W "a } ’ 


(51) 


which we write as 
E' - E q + 


E = - 


l T-0 + “ S } 


N 

e, = i Lim r 

^ nT-oJ 


f < xv ^ 


(52) 


The ground state energy is then 


E - E o + ®b + + E o 


(53) 


We note that E q is not the ground state energy of the electron 
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gas at density N/fi since the chemical potential p is that appro- 
priate to the complete system, namely electrons and ions. But 
has the same form as that derived in section III, for we 

may expand G(k,q,co ) in a Laurent series. 

/-0 \) 

, n ^( n ) 


= I X * G n 

n=0 


(54) 


so that, using (48), 


< XV > = |- Y Y V(-p)X n+1 G^ n ' ) (k,k-p;cD > )e°v 

A p fc=l /-i ^ r+s V 


CD 0 


(55) 


£,k,v n=0 

and the expression for E^ now reads 
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_ _ Lim 1 v” 1 2 V T7/ \p(n+l)« , s cd..0 + 

^"T-os 1 iMe 1 V( -E )G (Js>JS-p;“v> e v > 


(56) 


N 


n=0 


jg»k, v 


which is of the same form as (29) and constitutes a more formal 
derivation of it. 

In order to calculate E^ we need explicit expressions for 
the quantities G^ n+ ^. Considering the lowest order term we 


note that G 


(1) 


can be calculated in terms of known electron 


gas quantities. We have 


G^(k,k-p;cD ) = G^°\k,a> )V(p) A (k,co )G^(k-p,co ) (57) 


r*s V 


p v' 


os ^ v 


( 0 ) 


which is shown graphically in Figure 1. Here G (k,cD ) is 

v 

the Green function of the interacting electron gas and A (k,co ) 

22 ~ 

is the zero frequency vertex function. The second order 


term in the band structure energy is then from (56), 
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Lim 1 'Y 
T-0 g L 

k,p,v 

n*/ r^> 

which upon transforming to a contour integral gives 



|V (p) | 2 A (k,o> )G^ ' (k,o> )G^ 0 ' (k-p,o> ) 

1 1 D V ^ V /*//%/ \j 


(58) 


£<2) = I y 1 V (p) 1 2 f A (k,co)G^ (k,0i)G^°' ) (k-p, co)dD , (59) 

D 4_i J J3 ~ ~ ~ ~ 

p,k C 

r-f 

where C is the contour of Figure 2. From the definition of the 

22 

zero frequency dielectric function of the electron gas we 
therefore have 


■ (2) = i 1 


-I I v (p)I w (p) [e(p,o;p) " *■] 


(60) 


with w(p) defined in (7) and u being the exact chemical potential. 

The higher order terms in the expansion of G are, on the 
other hand, not well known, and the analogues of A^Oc, cd) must 
be approximated. 

We illustrate our approximation by re-calculating (58). 

Using the spectral resolution of G^°^(jg,co), i.e. 


G^°^P,cd) = ^- J <±d’ 


A(p,cd' ) 


(D-CD 


(61) 


we have 

(2) = 1 Lim 1 y f ^2 

d - T-0 p L J 2tt 2tt 

P,k,v 


V(p) A (k,o) ) 

1 1 p ~ v CD -CD- 

£ v 1 


A(k,co- ) A(k-p,co„) 


“v-“2 


( 62 ) 


which, exploiting a further transformation of the v sum to a 
contour integral gives two pole contributions, 
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do), dco 


„(2) 1 Lim v P 1 2 X ,2 W1 

^ ■ S T— 0 1 J l V( E^I A(k >mi )A(k- E ,m 2 ) 

N _ 


p,k 


A p (k,a> ) 

•{ ^ — n(co 1 ) + 

l -0^2 1 o^-to^ 


Ap(k,co 2 ) 


n(co 2 ) ]■ 


( 63 ) 


where n(o>) = (e^ ^ + 1) \ Our first approximation is to 

make an undamped quasi-particle ansatz for the spectral function, 
i.e. 

A(p,o>) = 2tt6 (<x> - d (p) - E, (p)) (64) 

rO Q ^ 1 ^ 

with s (p) defined to be the real part of the self energy 
satisfying Dyson's equation (f(j>) = cf o (p) + 2-^ (j0 * dT (p ) ) - 
Then (63) becomes, 

1 v - £o^ - *i<&» 

- 1 |V(P)I l [<f G (k) + E^k) -i 0 (k- E ) - S^k-p)] 

P,k 

(65) 

Ap(k, £ (k-p) + E-, (k-p))0(y. - C? (k-p) - E, (k-p) ) 

[ d_(k) + E n (k) - £ (k-p) - E-. (k-p) ] J 

Our second approximation is to replace Ap(k,co) in these expressions 
-1 23 

by e (p , o , p ) . We then have 

, „ , , eo - £ (k)-s (k))-e(u-f (k-p)-s .(k-p)) 

- L inp; > e(p,o,u) t (f (k)- £ (k-p))+(E 1 (k)-E.(k-p)) 

IN - U ^ O ^ ^ 1 ^ I ^ ^ 


P,k 


Furthermore, we write the chemical potential as 


( 66 ) 


u = p + 
eg 


9 


(67) 
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where yi is the chemical potential of an electron gas of 
density N/ft, i.e. 


U 

with k 


eg 

3 

f 


<~£ ^(^jd 

2_ 

3tt N/fi and 



( 68 ) 

h 2 2 

k^ so that (66) becomes 


r I W(E)I 2 T7?7^T{ e(£ f +6 u- - < s i< k ’ £< k » - 2 i< k f>^ g >» 

Vi t ^ 


p,k 


-0(<5 f +6^-<5 (k-p) - (S (k-p, d(k-p) ) - E, (k ,p ))) j 
x p i ^ ~ ~ it eg J 


(<f (k) - £ (k-p)) + (S,(k) - E,(k-p)) • 

o ^ o ^ ^ i ^ r-» 


(69) 


The final approximation is to neglect differences in self 

energies. For an electron gas at metallic densities this 

24 

approximation is fairly well satisfied. Thus the final 
approximate expression is 


- L e(p,o,p) { £00 " ( k ”P ) J 

IN «, r^f O'**' O r*t 


E>~ 


(70) 


For higher order terms we proceed in the same manner. Denoting 

the above approximation to V(p)A by a double broken line and 

~ P 

r— ' 

by a double wavy line the analogous approximation for the 
electron-electron interaction, we include the class of diagrams 
given in Figures 3 and 4. It can be shown that these corres- 
pond to a random phase approximation in the sense described by 

25 

Cohen and Ehrenreich provided one takes e(p,o,p) to be the 
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Lindhard dielectric function. 

We next examine certain complications which appear in 
fourth and higher orders and which are illustrated by the 
fourth order diagram of Figure 5. This gives a contribution 
to the band structure energy 

I 

fi I V(-k)A_k(E+k ) m v ) G(o) ( E>mv )v( a i)/V 3 i (p,m v ) 

JS>E>3l>3 2 ’ v 

• G ( ° 5 (P+9! ) V (g 2 ) A gz (p+ 3l .^ v ) G ( ° } (p+ 3i + 32 , o> v ) 

v (Js-ar32 )A k- 3l - 32 ( E + 3i + 32* m v )G(0)( E + !S> <B v ) • (71) 


In evaluating the v sum, we perform a contour integration and 
the possibility of double poles is evident (see Figure 6). 

The double pole contribution gives rise from differentiation 
of the factor [e^ ^+1] \ to a 6-function contribution in 


the T = 0 limit, i.e. 


AE 


-I I 

k,p,q 

r*/ 


V (k) 2 V (q) 2 

I I rO I 


G(k) 


I i 




6 ^f +6 u- W 


(£ n (p)“ £ n (p+k))(£ n (p)- £ n (p+q)) 


(72) 


From (B-5) and (B-8) the origin of this term is clear. It 

arises from an expansion of 9(E £ - E(k)) where E(k) is the 

r ~ ~ 

eigenvalue of the single electron band structure Hamiltonian. 

It is important to note that this expansion is invalid when k 
is near a zone plane: In fact AE of (72) diverges quadratically 
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there. Although the behavior of these anomalous contributions 

is general, we can ignore them provided the 9 functions 

occurring in the other expressions are modified from 9 (p^- <f^(k)) 

to 9(p^ - E(k)) where p^ = <f + 6^ = E r and is the chemical 
~ f p f 
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potential one computes in a band structure calculation from 

N = 2 £ 9 (E f - E(k) ) . (73) 

k 


The contributions from (71) not involving 6 -functions may 
be shown to give the first three terms of (B-8). The first 
term of this expression is well-defined, however, the second 
and third, due to the squared denominator, are divergent when 
the Fermi sphere is near a zone plane. This divergence is 
an artifact of the asymptotic nature of the expansion (B-l). 

In Appendix C we show that a resummation of diagrams leads to 
a finite result. 


Finally we make a remark concerning the electron gas term 
(p. ) - This can be calculated from approximate expressions 
for H(p) (e.g. the Nozieres-Pines formula). However, to gain 
some physical insight we expand H (p)+ pN about p = p , 

S 2 S 

3 0 (u) + UN = NE 0 (U 0 ) + i ( — -J ) (6|i b ) 2 + (74) 

3d 

and noting 



( 75 ) 
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we see that the right hand side of (74) becomes 


NW ' 2^T ( AN6 ^ b)2 + ••• (76) 

so that the change in electron gas energy lowers the total 
energy and is clearly related to the distortion of the 
spherical Fermi surface of the electron gas into the lattice 
symmetric Fermi surface of the periodic system. We may 
also observe that if (50) is written as 

S(u) = H 0 (n) + ^(li) (77) 

and expanded to fourth order in the external potential, the 
following expressions are obtained for internal energy, 
chemical potential and pressure: 


E ■ { W + + + } 

(78) 

+ 4 jf- K (6 u 2 ) 2 + 0(V 5 ), 

O 


b = 


where 


6p 2 

^3 

% 



(79) 

(80) 
(81) 
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, 1_ v . 1/2 d H 

+ K t 6u 2 9 V r s .2 


2 .(2) dE (2) 

- 2r A.) 


dr 

s 

r dB, 


s dr 


+ Q K T l 1 + 3 dr ), 

r\ o ^ 


( 82 ) 


and 


p “ { p o (p o ) + p(2)+ p<3)+ p(4) } + K t 6p 2 K r 

o L 


‘\ (2) 

s , 2 
dr 

s 


- 2r 


dE b 


( 2 ) 


s dr 


4 -V K, 

fi ■ 


,(6d 2 ) 2 (l 


1 

3 


r dB r 
s ] 

B dr 


0(V 5 ) 


(83) 


where 

(m) 

P = 


JL 

4nr 2 

s 



dr 

s 



(84) 


The quantity is the isothermal compressibility of 

the interacting electron gas and as before p o is the chemical 
potential of the interacting electron gas, both evaluated at 
density (The bracketed terms are to be expected from 

zero temperature perturbation theory.) We note that the two 
methods agree to third order but in fourth order differ for 
the physical reason outlined above (i.e., fermi surface dis- 
tortion). These differences although small are not negligible 
as will be shown in the next section. 

Recapitulating to this point, we have seen that the theory 
presented in section III must be modified in several ways. 
First, the electron gas term in the total energy must be 
corrected to take into account the shift in chemical potential 
due to the ions. Secondly, the expressions of section II for 
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except for the first, must be multiplied by an additional 
factor of e ^(k,o,p). Thirdly, terms such as 4b of Figure 4 
must be included in a self-consistent calculation. (These 
are essentially Hubbard's H diagrams which from his point of 
view are connected with double counting.) We now turn to a 
discussion of the magnitude of these various corrections 
for the particular case of a solid composed of protons 
arrayed on a Bravais crystal lattice. 


V. ATOMIC HYDROGEN 


In this section we present the results of calculations for 
zero temperature thermodynamic properties of three atomic 
hydrogen lattices, simple cubic (SC), face centered cubic (FCC), 
and body centered cubic (BCC) . This choice was made partly 
for convenience of computation, but more importantly because 
of the relatively large difference in Madelung constant 
between SC and the other two structures. We shall use the 
expressions (78) - (83) and proceed order by order. 

A. Electron Gas 

\ 29 

We have taken the Nozieres-Pines interpolation formula 

for the ground state energy of the interacting electron gas. 

,.q„ .2/3 , .1/3 .1 

W ■ Kf ) + (-0-115+0.031 fcr.). 

(85) 

In comparison of structures, the magnitude of the structure 
independent contribution plays no role so that a better approxi- 
mation is not necessary. In any case, the Nozieres-Pines' 
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expression compares very well with more recent forms 

B. Madelung Energy 

The Madelung energy may be written in the form 

®M = -V r s (86) 

where the Madelung constant A^ for the three structures is 
given by 31 : SC: 1.760122, FCC: 1.791749, and BCC: 1.791861. 
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C. Second Order Band Structure Energy 

\ 

We take the Lindhard expression for the dielectric function 
in the calculation of the terms in the band structure energy: 


*(”!%) = 1 r S ® (T|) 


1/3 


(87) 


g(Tl) 


= — rim 




■n 


2 L 4ri ' rr 


with ri = k/2k£. 

Then the second order band structure energy may be written 

.( 1 , 2 ) 


as 


E b (2) <%> " - 72 I 

6,1 TffO 


_glnl 


-2 

11 ,,174 \ 1/3 t . 

1 + \% ) r s 8 (Ti) 


( 88 ) 


D. Third Order Band Structure Energy 

I 

This contribution is given by (B-7) and corresponds 
to the diagram 3 of Figure 4. It may be written in the 
following form: 


E b 3)(M o ) = ■ r s Xw(-T J )w( 3 -T 1 i)«(rji)H (3) (T),t J i) (89) 


16/ 4 


1/3 


(3) 


■n, ri 


where w(ri) = — ~ 

Z 


, and 


i\ e(^l,u ) 

O 


H^On.np = k £ | 


n „ (q) 


' ^o ( 3 + !S )J[£ o ( 3> " < ~(3 + ~pl 


(90) 


,(3) 


The complete expression for H v y (r\, ) is given in appendix E. 
We note here the following asymptotic properties: 
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u(3) / „ X 1_ 

H (TIjTl,) o 

r* ~ «, , Z. 


16n‘ 


Till- 


->/w 


1 -1/ 1-a 

tan ./ Tta 


A /\ 

a s TV TV, 


> ri,^ « 1 


(91) 


» (3) - rh -T2 i 1 + 1 (-T + -T + ^7 a )} • 


(92) 


96tt r\ p 


TV, 


ri,'n 1 » l 


The third order contribution thus depends linearly on r apart 

s 

from a weak dependence contained in the dielectric functions. 
(3) 

The function H v J (11511-) in this approximation is independent 

rsj rvi 

of r and depends purely on the structure. It is everywhere 
s 

finite but has discontinuous derivatives for certain values 

of T|, as discussed by Lloyd and Sholl^. We have expanded 

(3) 1/4 \ 

(p Q ) as a power series in the parameter cr s = 9^ ) r s 

which occurs in the Lindhard function, 

E^ 3 ) (u 0 ) - ar s (b D + cr s b x + (cr^ 2 ^ + (cr g ) 3 b 3 + ...)(93) 


1/3 


where a = - 9^9^" ) • The values of these structural 

constants are given in Table I. 

E. Fourth Order Band Structure Energy 

There are several distinct contributions in this order. 
We consider first the most divergent parts of the last two 
terms in (B-8) namely 

2 L n o^ L 177k 

'i 


2 v - V(K ) 

-|XVS> 

k aio 


<<s 0 ~ £ t? 


(94) 
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and 

- il 

k 

r^> 

which we write as 

E 4'I[ E 1 4) <V + E 2 4)( V] 

i 

The explicit expressions for and E 2 are given by 



V(K ) 

~J6 


6 < E f£o> 


L £(K,)' ,/ c .2 


( 95 ) 


(96) 


(97) 


E 


(4) 


w-i^(£) 2/3 4^] 4 i 


i° „ 2 . 

*n (l-n ) 


(98) 


In Figure 8 we show E-j^ (r\) /E^ 0 ^ (ti) and E^^ (ri) /E^^ (tO as 

functions of T] (where = S'E^^'n)) along with the resummed 

“9 (4) 

expression given appendix D. Note that E^ On) is part of 
the anomalous contribution as discussed in section IV and that 
it must be included at finite order to give the appropriate 
limiting agreement with the resummed diagrams. Furthermore 
we note from the positions of the first reciprocal lattice 
vectors that the contribution of this term will be small. 

The behavior exhibited in this term is representative of the 
nature of any spurious divergences introduced by zone planes 


and illustrates the interconnection between band structure 
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effects and the methods (finite T and T = 0) of perturbation 
theory. 


Secondly, we consider contributions from the diagram 4b 


of Figure 4. This term may be written 

(4b) _64 /_4_ \ 3 y' 1 g 2, . 

E 2 27tt \9tt / r s L 2 , , G ( V 

V (tip 


(99) 


where 


i 

G ^ T 'l^ = X ~2 


„ r '2 €(t ' 2 ) ( 3l'2 2 ) e( 3l’3 2 > 


{ 2h(3) (Z!i>3 2 ) +h3 (3 2 ’I! 2 -9x>] 


( 100 ) 


and can be calculated readily since the expressions for H' '(p , ri 0 ) 

'-'Z 

are known. Furthermore, apart from the weak r g dependence of £ 

3 

this term is proportional to r g . Numerical values for two 
representative values of r g are given in Table II. 

We next consider the correction which arises as a 


consequence of the chemical potential shift, namely the last 
term in (78). 

E i 4b) = $ ir k t (6m 2 )2 ’ (101) 

o 

This is known from the expressions for the compressibility 


of the electron gas and the second order value of the chemical 
potential. In fact, as a consequence of the compressibility 
sum rule, it may be shown that this term is precisely given by 
the diagram for in the limit that the momentum transferred 

by the internal Coulomb line approaches zero. 



Finally we consider contributions due to diagrams of the 


form labelled 4a in Figure 4. There are two contributions 

apart from those already discussed in the first part of this 

section and are given in (B-8). One is an off diagonal part 

V (K ) V(-K.) V(K.-K) i 

~x Vi 1 

;J e (K. ) e (K. -K 


n (k) 
o ~ 


o 0 v ' ; ^.i iv 

° i *(K ( ) e(K.) e(K,-K ( ) * 0 L (^ V (< V 9 (<f o" 9 


i#),A 

j^O,A 

J#0 


(102) 


and the other has diagonal parts: 

V (K ) 2 V (K . ) 2 -i 

2Q o I ^7(KjTi I n o ( y{ ^ ..2 

tUn „ 1, (<- " 


i^O, a 

jl¥ 0 


o w jt' ^o”^i-P 


and 


V (K ) 2 V (K. ) 9 6(E S 6J 


I V 1 -i 1 V v 

"N I e (K ) ‘ e (K . ) L 

• un „ ~X i. O "X O X 


i^O, A 

liO 


< W 


} 


(103) 


The second part of (103) is an anomalous contribution, which 
disappears along with the singularities from the double poles 
if the resummation of appendix D is used. These terms are 
awkward to handle in numerical work, although in principle 
there is no difficulty. [One problem is the time needed to 
calculate a nine dimensional sum. Another is that the kernel 


K(p ,ru,t) 


>4 1 


k 


n 00 

o ~ 


(f.-i.xvy'vv 


(104) 


has, as yet, no analytic representation. We have been able to 



reduce it to a two dimensional integral. It has asymptotic 
expansions which give 


k (tJi>t]2» 



1_ 1 

48tt 2 (2k f ) 3 


111, , 

2 2 2 r >2 ,r '3 ^ 1*^) 

ni ti 2 t, 3 


V 2 


1 ^/ A a A \ 
f ri 2 » ^3) - 


1 r-/ A A A N 1 

— £(^.113.^) - — 


2 ’3 


VV^ « 1 . 


f (^2,113,11-,^) 


(105) 


where f is a complicated function of angles. We calculated 
these terms (102 and 103) by taking as an approximation for 
K( t 1i > r lo>'n 0 ) its large p expansion and by setting /- t = 1. 

'-'■'1 r^i. G fH/ 

The former is an underestimate but note that for the structures 
we consider r\ is always > 1. The latter is an overestimate. 
The form is then: 

.(4a) ~ 1 ft\ Z/3 ,2f JJVj: 3 , 


4(3rr) 


2 \9tt 


) 'si 


•Hi 


2\L 4 , v 2. 
\ 2 ^2 ^1"^) 


* ■ 


r c, 
s 4 


31 ~2‘ 


(106) 

2 

which is proportional to r and is probably an underestimate 

s 

overall. The values for the factor c^ are given in Table I.]. 

In Tables III-V, we give the thermodynamic functions p,E,G, 

at T = 0°K calculated to third order in the electron-ion 

interaction. In Table II, we list the explicit contributions 

to fourth order at r =1.6 and r =1.36 corresponding to 

s s 

low pressure and 1.9 Mbar respectively. 3 ^ The contribution 
is an estimate as noted above. Note the approximate cancellation 
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in the fourth order, and further that at high pressures the 
SC lattice is predicted to be unstable relative to FCC and 


BCC. 



VI. DISCUSSION AND CONCLUSIONS 


We have given a procedure for calculating the ground 
state energy of a simple metal and have shown that there are 
basically four contributions involved, viz. electron gas, 
static dielectric energy, Madelung, and core exclusion. 
Furthermore, we have seen that the shift in chemical potential 
from that of a uniform electron gas must be taken into account 
in calculations going beyond 2nd order. In particular, we 
have emphasized that T = 0 time dependent perturbation theory 
does not give the true ground state when the unperturbed 
system is taken to have a spherical Fermi surface (a fact 
first noted by Kohn and Luttinger"^) and have shown the relation 
ship of this to the deformation of the unperturbed Fermi 
surface. We have observed that if one expands the free energy 
uniformly in powers of electron-ion interaction, differences 
between finite and zero temperature perturbation theory appear 
only in fourth and higher orders, and furthermore that certain 
divergences at zone planes can be resolved by resummations. 

The preliminary calculations reported here for atomic 
hydrogen seem to indicate that a happy cancellation may occur 
in the fourth order, at least for theSC, FCC, and BCC structures 
although more detailed calculations are required to be certain 
of this. The calculations reported have been done using the 
Lindhard dielectric function. In third and higher orders this 
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is a very good approximation since the dielectric function 
occurs as 1/e. However, in the second order e - 1 appears. 

This acts to change the magnitude of the second order contri- 
bution slightly but does not affect the energy differences 
between SC and the two close-packed structures. The use of 
the Lindhard function, as noted in section III corresponds 
to a self consistent Hartree (RPA) approximation. We remark 
that the zero pressure density of the structures studied will 
be extremely sensitive to the exact fourth order corrections 
due to the weakness of the minimum in the free energy as seen 
in Table IV. Also, a third order calculation predicts an 
instability of the SC structure relative to the two close 
packed lattices at a pressure of ~2 - 3 Mbar. The 
exact transition pressure is again sensitive to the magnitude 
of the fourth order corrections. It is clear, however, that 
such a transition must appear at some pressure, for the band 

structure corrections depend upon positive powers of r whereas 

s 

the Madelung term depends inversely upon r . Thus eventually, 

s 

the static lattice having the lowest Madelung energy should 
be most stable. 

Brovman et al.^ have computed ground state energies for 
atomic hydrogen at zero pressure by using the T = 0 expansion 
to third order in the electron-ion interaction and found an 
interesting class of low energy anis tropic structures. We 
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regard the effect of fourth order corrections to these calcu- 
lations as an open question, but one that can be settled using 
the above expressions. Finally, we again emphasize that we have 
treated the lattice as static and that it will be necessary 
to consider lattice zero point energy in a complete determination 

of structural stability since the zero point energy is of the 

( 3 ) 

magnitude of E£ . Calculations of such phonon effects are 

33 

in progress. 



Acknowledgments 


We would like to thank Drs. B. Nickel and A. 
for numerous helpful discussions. 


- 39 -. 

I 


B. Bringer 



Appendix A. 

We derive (29) by writing W in the form 
? = Y K (k) V(-k) V (k) + Y K (k,q)V(k + q)V(-q)V(-k) 

^ l L » i r** ^ rw L i Z. rv/ />-» /%/ 

is ~’3 

+ l K 3 ( iS>3l>S2> V( 3 + Sl + a2 )V <-y V <-S 1 ) V <-32 ) + 
S*si,a 2 

Defining symmetrical coefficients by 
L x (k) = [ K x (k) + K 1 (-U)l 

L 2 (k,q) = [K 2 (k,q) + K 2 (-k-q,^) + K^C^k)] 

= t K 3 (k,qi,q 2 ) + K^-k-q^q^q^q^ 

+ + V&’Sry 3 


(Al) 


(A2) 


• > 


we may write (Al) as 

? = k Y L (k)V(k)V(-k) + i Y L.(k,q)V(k+q)V(-q)V(-k) 

Z L-i nss /*-» O L-i Z r* 

k k,q 

r>* rw* />-» 

+ II L 3 ( ~’3i’S 2 ) v( iS + 3i + 3 2 ) v(_ y v( '3i ) v( '32 ) + 

k,qx»92 

Taking the variation of (A3) with respect to V(-k) we have 

5 1 H(-k) - 1 L i<is> v <y + 1 l 2 <M> v <!s + s> v <-s> 

k ~ k k,q 


(A3) 
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+ £ L 3 (k,q 1 ,q 2 ) V(k+q 1 +q 2 ) V(-q x ) V(-q 2 > + 

!S»3i»S2 

Equating (26) and (A4) we have, noting that the V's are 
arbitrary external potentials, 



and hence (29) • 


(A4) 



Appendix B . 


To derive an expansion of (41) we write the band energy 


as 




E(k) = C (k) + Y 

O ^ L-X 


V .V. 
ox 10 

if - £. 
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(Bl) 


and the occupation number as 

n(k) = n (k) - 6 (E £ (k))- Y 

O ^ II O ~ 


V .V. 

Ol IQ 

r V4 
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(B2) 


where n (k) = 9(E £ - £ (k)) and n(k) = 9(E,--E(k)). Clearly when 

0 ~ t o ~ ~ r ~ 

k is near a zone plane these must be viewed as asymptotic. We 
find the following expressions for the Fourier components of the 
density 
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Using (29) and supplying the extra factor of e \k) in the 
third and higher orders, we find for the energy, 

e (2) - I V , V / K ,|2 1 
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Appendix C . 


The diagrams which correspond to the second and third 
terms of (B8) are shown in Figure 6. The two diagrams of 
6a are equal in magnitude when summed over k,(j so we need 
only calculate one and multiply the result by a factor of 2. 

We now observe that the series of Figure 7a may be summed, i.e. 


VE’V I [ 'w k l 2 G 0 (p+!s.» v )G 0 (p.» v )] 
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(P»^ v ) 
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Hence the series of Figure 7b can also be summed, and supplying 
the factor of 2, the resummation gives for Figure 6a, 


I I \\\ J x dx ^r 
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which no longer has double poles and hence is always finite. 
Similarly, the contribution of Figure 6b is the first term 
in the series of Figure 7c which may be summed to give, 

N I \\\ J d(X ){“7T -1 2,~ .2 
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This again has only simple poles, and moreover is seen to be 
a correction to rather E^^ . In fact, the integrals 

appearing in (C3) can be done analytically. 



Appendix D . 

We give an explicit expression for (C3) by noting that 
the roots of the equation 

G“ 1 (p,co)G" 1 (p+k,co)-X 2 |V | 2 = 0 (Dl) 

O r*-> o K. 

r+s 

are given by 

£ + - fc{(%)+<f( B -k)) ± t£<p-kW(p)) 2 + 4x 2 |V k | 2 ]*}, (D2) 

where £(p) = & (p) + S, (p) as in (64). Then (C3) becomes 

r-> O r- J_ 




^ n — 


k 

with 


0) (jo -a> ^(p)-4 (p-k), 

*T* “ T u O ~ O ~ ~ 




(D3) 


“+ = i{ (? o ( E )+,f o ( E-y> ± f( ,f o ( E-y- £ o ( E ))2+4x2 l 9 kl 2ji } 

r-~> 


(D4) 


and where we have assumed the self energy to be momentum 
independent as discussed in section III. We rewrite (D3) as 


I rvy - E 4 0 ) <y ] 


(D5) 


where E, (k) is given by the first two terms and E^ 0 (k) by the 
4 ^ 4 ^ 

last. Then in terms of the dimensionless variables: 


% 2k, 


v = 


Vk 



(D : 


L 


= _a. 


(D 


(D6) 


A = ^[44 + v ] ? 


The integrals may be evaluated and we find for E, (tc) 
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(i) 8 < i - v 


vy - 


3/ 9rr_ \ 2/3 

" 4\ 4 ) 2 

r 

s 


j^4(3£-4+A ')—j~ Qm\ 


[6+1+2A] g +A+l | 
[d+l-2A]^-A+l 


+ [£ + 1+2A]^[-^|^- [ 3£-4+A 2 *-2 (1+A) 2 ] + |(£+1+2A) 


+ ^ (£.+l+2A)(£+l-3A) + 4(A 2 -£)] 

- [t+l-2b]^j^- [3£-4+A 2 -2(l-A) 2 ] +| (£+l-2A) 

4 o i 8 5/2 

+ ^ (S+1-2A) (£+1+3 A) + 4(A -6)] - ~C } , (D7) 


(2) C> 1 - v; £< 1 + v 

E 4« - - Kf - ) 2/3 "T {i(3^-4+A 2 ) ^ fa| . [^ 1 ± 2 AjV 1 | 

4 4 r 4 [A -£] 2 

s 

+ [e+l+2A]^[-^^-(3£-4+A 2 -2(l+A) 2 ) + ^ (6 +1+2 A) 

4 2 8 5/2 , 

+ g (£+l+2A)(t+l-3A) + 4(A -£)] - £ }(D8) 


(3) £ > 1 + 


v 


e 4« - - Kt) 2/3 ^{i(3f-4+A 2 ) 

4 [£+1-2 A] 2 +A-l 

b 

+[£+l+2A]^[-^^- [3£-4+A 2 -2 (A+l) 2 ] + |(£+1+2A) 


+ ^(£+l+2A)e+l-3A) + 4(A 2 -<0] 
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- [6+l-2A] 2 [-^^-[3£-4+A 2 -2(A-l) 2 ] - | (5+1-2A) 


(6+1- 2 A) (6+1+3A) - 4(A Z -e)] - } (D9) 


5/2 


15 


and for (k) 

4 ~ 


E l 0) <y - - b(t ) 2/3 -T v 2 {(£-l)fcl^Tll + 2/£ } . (DIO) 


To see the relative magnitude of these contributions we 
compare them with the second order contribution which we write 


as 


(Dll) 


^ 2) 'I E 2« ' 

k 

We distinguish between the cases in which E 2 is evaluated 
using e(k,p ) and e(k,p) by E^°^ and E 0 , where the u dependent 
Lindhard function e(k,p) is given by 

£( ~ u) “ 1 + , 9n s .i/3 ^ + x ] (D12) 


Hf) 


with 


T1 = 


b 2k j 




2/p 


b * 


In Figure 8 we show [E^ (ri) -E^°^ (ri) ] /E^ 0 ^ ("n) with p = p^. 



Appendix E . 


The function (90) was originally evaluated by Lloyd and 

g 

Sholl . We write it in a somewhat different form: 


h<3) <VV ‘ :^2 


P 2 +1 


- .27 Re {< r| 1 - ,1 2 COS0 )lo s (ri7^T ) 

64rr r^r^sm 0 '2 

■n-,+1 


(El) 


+ (p 2 - T ii cos e)i°g(^TT ) + 3 iog[ (r^-i) (ri 2 -i)] 


r 2 2 2 2 i 

- p log|[n 1 'n 2 +ri 1 +-n 2 - A-ri^r^cosG + cos20] +^(^ 112-008 0)1 


where 


2 2 2 - 

P = (■n^+'n 2 -2Ti^ri 2 cos0-sin 0) 2 

When p = ip 1 , i.e. when ru , -T| 0 form a triangle inscribed 

'"'1 

in a circle of diameter 1 this function becomes 

r\ 2 +i n 1 +i 

|(il 1 -ri 2 cos0)log|— ■ j\ + ('n 2 -p 1 cos0)log|— 


2 2 
64tt ry^-r^sin 0 


2 2 2 2 l l 

+ p 1 arg[n 1 i 2 +ri 1 +'n 2 -4ri^ri2 cos O +cos 20] + 2p ' ('n 1 'n 2 -cos0)| j 

with p’ = (sin 0 - t^-t^ + 2p^ri 2 cos0) 2 and the argument function 
is the principal branch with the branch cut along the positive 
real axis. 
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32. For comparison we have calculated the energy to second 


order for the liquid metal from 



where S is the hard sphere Percus-Yevick structure factor 
JLi 

and e u is the Hubbard dielectric function modified to 

n 

satisfy the compressibility sum rule. For r =1.6 

s 

and r\ = packing fraction = .45 we find E = -.976, these 
results bong very insensitive to r\ (~2 % for t)-''. 7 - .3). 

Most of this change can be traced to the Madelung energy. 

33. A correct evaluation of the phonon spectrum is essential 
in determining the zero pressure density for a given 
structure, as well as in assessing dynamic stability. 
Futhermore, only when the presence of phonons is taken into 
proper account can the virial theorem (E = -K + 3pfi, 
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where K is the kinetic energy of electrons and ions) 
be satisfied. 



Table Captions 


Table I : 


Table II: 


Table III: 


Table IV: 


Table V: 


Parameters in expansion (93) of third order band 
structure energy. 

Contributions to fourth order in electron-ion 
interaction to free energy. - Fig. 4b, 

E^ a - Fig. 4a, E^ - Fig. 7c, E^ b ^ - chemical 
potential correction - see text of Section IV. 

T = 0°K equation of state for atomic hydrogen 
(to third order). Note that these results are 
appropriate to a static lattice and do not, there- 
fore, include phonon contributions to the equation 
of state. 

Free energy at T = 0°K for atomic hydrogen vs. r 

s 

(to third order). 

Gibbs free energy at T = 0°K vs. pressure for 
atomic hydrogen (to third order). 
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Table I 


Lattice (real space) 

b 

0 

b i 

b 2 

b 3 

C 4 

SC 

.08202 

.1195 

.1506 

.1748 

-.00310 

BCC 

.06483 

.06591 

.05467 

.04050 

-.00275 

FCC 

.06663 

.06945 

.05933 

.04555 

-.00260 


- 56 - 














-57 




Table III 



NOTE : 1 unit of pressure = 147.15 Mbar 


- 58 - 















































Table IV 


r 

s 

free energy 


SC 

FCC 

BCC 

1.65 

-1.04803 

-1.04338 

-1.04209 

1.64 

-1.04807 

-1.04353 

-1.04224 

1.63 

-1.04805 

-1.04361 

-1.04233 

1.62 

-1.04796 

-1.04363 

-1.04236 

1.61 

-1.04781 

-1.04360 

-1.04233 

1.60 

-1.04759 

-1.04345 

-1.04222 

1.55 

-1.04538 

-1.04188 

-1.04062 

1.50 

! 

-1.04104 

-1.03818 

-1.03693 

m 

-1.03414 

-1.03197 

-1.03073 

mm 

-1.02414 

-1.02272 

-1.02149 

1.35 

-1.01042 

-1.00979 

-1.00858 

1.30 

-0.99217 

-0.99242 

-0.99122 | 

i 

1.25 

-0.96842 j 

-0.96961 

-0.96843 | 

\ 

1.20 

-0.93796 

-0.94019 | 

-0.93902 j 

1 

1.15 

' 

-0.89928 i 

-0.90262 

-0.90147 j 

t 

1.10 

-0.85045 

-0.85502 

■ 

-0.85388 

1.05 

j 

-0.78903 i 

-0.79495 

-0. 79383 

o 

o 

• 

t—i 

-0.71188 j 

1 

-0.71929 ; 

-0. 71819 
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■ Table V 


pressure 

Gibbs free energy 


SC 

FCC 

BCC 

0.0 

-1.0481 

-1.0436 

-1.0424 

■HE 

-1.0390 

-1.0349 

-1.0336 1 

1.0 X 10" 3 

i 

-1.0305 1 -1.0266 

-1.0253 

5. 

-0.9707 

-0.9683 

-0.9670 ! 

1.0 x 10“ 2 

i 

-0.9092 1 -0.9080 

-0.9068 

2. 

-0.8081 : -0.8 085 

-0.8073 

3. 

i 

-0.7233 i -0.7248 

-0. 7237 

5. 

-0.5809 

-0.5841 

-0.5829 

1.0 x 10' 1 

-0.3019 ! -0.3085 

-0.3075 

5. 

1 

-0.9707 ■ -0.9683 

-0.9670 

1.0 

i 

1.8572 • 1.8377 

1.8387 

5.0 

# 

5.6614 i 5.6273 

5.6282 
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Figure Captions 


Figure 1: 


Figure 2: 
Figure 3 : 


Figure 4: 
Figure 5: 

Figure 6 : 
Figure 7: 

Figure 8: 


Figure 9 : 


First order correction to the Green Function. The 
solid line represents the electron gas Green function, 
the dashed line the bare external potential, and 
the triangle the vertex function of the electron gas. 
Integration contour for (59). 

Corrections to the Green function. The double 
dashed line and double wavy line represent the 
dielectric approximation to the vertex function 
described in the text. 

Contributions to the band structure energy. 

Fourth order contribution to the band structure 
energy given by (71). 

Divergent fourth order diagrams. 

(a) Partial summation for Green function. 

(b) Partial summation for the diagrams of 6a. 

(c) Partial summation for the diagrams of 6b. 

[E 4 (p)-E^° ) (p)]/E 2 (o) (r 1 ) (cf. Appendix D) . 

— ^(ri)/E^ o) (ti) 

[E{ 4 ) (p)+E^ 4 ) (r 1 )]/E^ o) (p) (cf. (97) 

and (98) and Appendix D). 

Note that left hand axis is right hand is p. 

Gibbs free energy difference relative to the simple 
cubic lattice for FCC and BCC metallic hydrogen. 
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Order 
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Fig. 4 
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Fig. 5 
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Fig. 7 
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